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Abstract 

In some quantum gravity theories, a foamy structure of space-time may lead to Lorentz invariance violation 
(LIV). As the most energetic explosions in the Universe, gamma-ray bursts (GRBs) provide an effect way 
to probe quantum gravity effects. In this paper, we use the continuous spectra of 20 short GRBs detected 
by the Swift satellite to give a conservative lower limit of quantum gravity energy scale Mqq. Due to the 
LIV effect, photons with different energy have different velocities. This will lead to the delayed arrival of 
high energy photons relative to low energy ones. Based on the fact that the LIV-induced time delay cannot 
be longer than the duration of a GRB, we present the most conservative estimate of the quantum gravity 
energy scales from 20 short GRBs. The strictest constraint, Mqg > 5.05 x 10^^ GeV in the linearly corrected 
case, is from GRB 140622A. Our constraint on Mqg, although not as tight as previous results, is the safest 
and most reliable so far. 
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I. INTRODUCTION 


A foamy^structure of space-time on short time and distance scales has been proposed in quantum 
gravity [ll-y]. The non-trivial space-time can lead to the violation of Lorentz invariance at the 
Planck scale. Probing the Lorentz invariance violation (LIV) effects provides a useful way to test 


the validity of quantum gravity theories, such as loop quantum gravity [ 4 , 5|, string theory 
double special relativity Q]. According to the LIV effect, massless particles have energy-dependent 
velocities. Hence, the velocity of a photon propagating in vacuum may have a tiny deviation 
from the trivial value, c. The delay of arrival time induced by the LIV effect is a monotonically 
increasing function of the photon energy and the distance of source. We need plenty of distant 
energetic photons in order to observe the LIV effect. 

As the most energetic explosions in the Universe, gamma-ray bursts (GRBs) can be detectable 
up to distances as far as tens of Gpc away from us. It has been proposed that GRBs provide an 
effective way to probe the LIV effect because of their cosmological distance and rapid emissions of 
energetic photons 0, 1^. In fact, GRBs have already been widely used to constrain the LIV effect 
[ 9 - 29 I. However, since we have little knowledge of the emission mechanism of GRBs, we cannot 
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and 


distinguish the LIV-induced time delay from the intrinsic time delay. Ellis et al. m performed 
a linear regression analysis of GRBs with measured redshift. The data were htted by a straight 
line with a slope corresponding to the quantum gravity scale and the intercept representing the 
possible intrinsic time delay inherited from the sources. They found a strong correlation between 
the parameters characterizing an intrinsic time delay and a distance-dependent propagation effect. 
Their work based on the assumption that all GRBs had the same intrinsic emission mechanism 
and the intrinsic time delay. Since the durations of GRBs span about 6 orders of magnitude, there 
is no persuasive reason to believe that two high energy photons from two different GRBs (or two 
photons with different energy from the same GRB) have the same intrinsic time delay relative to 


the trigger time of low energy photons. As an improvement, Zhang and Ma 


211 fitted the data 


with straight lines of the same slope but different intercepts. However, photons on the same line 
must have the same intrinsic time delay. This is not always true, because photons on the same 


line may come from different GRBs and have very different energy. Chang et al. 


Ill ] used the 


magnetic jet model to estimate the intrinsic time delay between emissions of low and high energy 
photons. Unfortunately, the magnetic jet model depends on some unobservable parameters, and 
thus introduces large uncertainties. Vasileiou et al. [l3 adopted maximum-likelihood analysis to 
test quantum gravity model in which photons had normally distributed velocities. They used the 
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synthetic data of GRB 090510, and chose the time interval and threshold energy ilith according to 
increased sensitivity and minimal systematic biases. However, in their calculation, only hundreds 
of photons below E'th and tens above were used. 

The previous works can be divided into two major classes. The first class constrains Mqg 


using one or some isolated high energy photons from GRBs 


Ill4l3l|. The second c 


Mqg using the observed spectral lag between the low and high energy band 


ass constrains 


IJ-lSQ]. The first 


method is not statistically significant because the number of high energy photons is too small. 
More importantly, we cannot know if high energy photons and low energy ones are emitted si¬ 
multaneously. The second method is based on the fact that, in some Fermi GRBs, high energy 
photons show a systematic spectral lag with respect to low energy photons. By assuming that the 
LIV-induced time delay cannot be longer than the spectral lag, Mqg can be constrained to the 
order of Planck energy scale. This method has the underling assumption that the low and high 
energy photons are emitted in the same region, which is widely accepted to be true since low and 
high energy photons form the same types of spectra (such as the Band function). However, there 
is still a possibility that low and high energy photons come from different regions. For example, 


Kumar & Duran 


30 1 showed that the high-energy photons (> 100 MeV) of GRB 080916G were 


unambiguously generated in the external shock via the synchrotron process, while the lower energy 


photons had a distinctly different source. Besides, Ellis et al. 


16l | analyzed the light curves of 


35 GRBs and found a systematic tendency for more energetic photons to arrive earlier than low 
energy photons. This conflicts with most Fermi GRBs which often show that high energy photons 

With the present knowledge, we cannot distinguish the 


come later than low energy ones 


31-; 


LIV-induced time delay and the intrinsic time delay from the observed spectral lag. If high energy 
photons are intrinsically emitted earlier than low energy photons (we can’t exclude this possibility 
even high energy photons are observed later), the LIV-induced time delay may be much longer 
than the observed spectral lag. On the other hand, most short GRBs show no significant spectral 
lag. If we still use the spectral lag as the upper limit of LIV-induced time delay, we may constrain 
Mqg to be infinity. But we cannot safely say that the LIV effect does not exist. 

In this paper, for the first time, we use the duration of short GRBs as the upper limit of the 
LIV-induced time delay. This is always true, regardless of whether the observed spectral lag is 
totally due to the intrinsic or LIV-induced time delay (or the combined contribution of these two). 
We use 20 short GRBs with measured redshift from Swift Q to give a conservative lower limit 
of quantum gravity energy scale. Whether the intrinsic emission process contributes to the time 
delay or not, the LIV-induced time delay is certainly smaller than the duration of a GRB. The 
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rest of this paper is arranged as follows. In Section II, we review the observational properties of 
GRBs. In Section III, quantum gravity energy scale is constrained using short GRBs. Finally, a 
short summary is given in Section IV. 


II. OBSERVATIONAL PROPERTIES OF GRBS 


Observations of GRBs, especially by means of the Compton, Swift, and Fermi satellites, con¬ 
tribute to the research of their properties. One of the most important properties of GRBs is the 
duration. The duration of a GRB is usually characterized by Tgo, during which from 5% to 95% of 
the total photon events in a specific energy band are detected. However, the duration Tgo depends 
on the sensitivity of the detector and the energy band in which the detector works. A detector 
with a lower and broader energy band generally gets a longer Tgo for the same GRB 351]. The 


observed durations span about 6 orders of magnitude, from milliseconds to thousands of seconds. 
The distribution of durations is bimodal and separated at around 2 s. Thus, Kouveliotou et al. 
3a | proposed a classification according to Tgo: long bursts with Tgo > 2 s and short bursts with 
Tgo < 2 s. The GRB light curves are irregular. Some are variable with many peaks, while some 
are smooth with simple temporal structures. 

The GRB spectra are non-thermal. The energy of the prompt emission is concentrated in the 
hundreds of keV range. In some of the brightest GRBs, photons with energy higher then 100 MeV 
(and maybe up to tens of GeV) have been observed [l^, [SlnSStj. X-ray emission is weak, and a 
small part of the emissions are below 10 keV [^. A typical GRB spectrum can be fitted with the 
so-called Band function 


38| |. It can be written as 


N{E) = < 


E 


100 keV 
100 keV 


E 

exp I - — 
Eg 


exp (/3 — a) 


E 


E < {a — 13) Eq 

A , E>ia-P)Eo 


( 1 ) 


100 keV^ 

where N{E)dE is the photon number in energy interval dE, a and /3 are the spectral indices, and 
To is the break energy. Typical evaluations for the spectral parameters are a ~ —0.5 ~ —1.5, 
^ ~ _2-3, To « 0.1 ~ 1 MeV 


39|, 


40(1 . Eq.([T]) is a time integrated spectrum of photon number 


and consists of two smoothly-joined power law parts. The Band function peaks at around a few 
hundred keV and it fits well to most of the observed spectra [35 j. 


The Swift satellite, launched in November 2004, has three instruments working together to 
observe GRBs in prompt emission and afterglow phases in the gamma-ray. X-ray, ultraviolet, 
and optical wavebands. The Burst Alert Telescope (BAT) first detects the GRB and accurately 
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determines its direction in the sky. Then in less than approximately 90 seconds, the X-ray Telescope 
(XRT) and UV-Optical Telescope (UVOT) slew to the GRB and start observing. Swift provides 
the light curves of prompt emission and afterglows and evaluates the duration Tgo of a GRB. In 
the following, we use photons detected by BAT, which works in the 15 — 150 keV energy band. 

III. CONSTRAINTS ON QUANTUM GRAVITY ENERGY SCALE 


In quantum gravity, there are several possible scenarios with respect to the breaking of Lorentz 
invariance. In consideration of experimental tests, Amelino-Camelia and Smolin ^ sorted them 


into three broad categories, i.e., naive Lorentz symmetry breaking, Lorentz symmetry breaking 
in effective field theory, and doubly special relativity. In all three models, Lorentz symmetry is 
broken at a very high energy scale, which is expected to be around the Planck energy. I n ge neral. 


the deformed dispersion relation of massless particles at leading term can be written as 


42| 




1 ±(^ 

\MqG 


( 2 ) 


where n = 1,2,3 denotes the linear, quadratic and cubic corrections to the dispersion relation, 
respectively, and Mqg represents the quantum gravity energy scale. The ‘+’ or in Eq. (i2]l 
corresponds to the superluminal or subluminal motion of particles. When the energy of a particle 
reaches a value at the same order of Mqg, the LIV effects will become obvious. We use the 
case where particles with high energies travel slower than those with low energies. 

As a result of the LIV effect, two photons with different energies, emitted simultaneously from 
the same GRB, have different travel speeds and arrive at the earth at different times. Due to the 
cosmological distance and the high energy of GRB photons, the time delay is accumulated as the 
propagation of photons and should be measurable. Taking into consideration the expansion of the 
Universe, we get the expression of the LIV-induced time delay of two GRB photons [^ . 43]: 

1 + n f AEY r , , {l + zT 


AtLIV = 


2i?Q V M( 


dz' 


(1 + z')^ + Qa 


( 3 ) 


where AE is the energy difference of two photons, and z is the redshift of the GRB source. Hq is the 
Hubble constant, and Qm and Da are the present values of the matter density and cosmological 


constant density, respectively. Throughout this paper, t 
Hq = 67.8 km s“^ Mpc“^, flm = 0.308, and Qa = 0.692 


re Planck 2015 results are used, i.e.. 


M\. 


As was mentioned in the introduction, some works constrained LIV using a small number of 
GeV photons selected from a few GRBs. The energy spectrum of a GRB at the GeV scale is 
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usually discrete. Only a few GeV photons can be detected in a GRB emission [4^. The results 
of using GeV photons to constrain LIV are unpersuasive in statistics. However, at the keV order, 
the energy spectrum is continuous, and thousands of photon events can be recorded. Thus, using 
keV photons from GRBs to constrain LIV is more statistically reliable. Taking GRB 090510, for 
example. Figure [D shows the light curves of GRB 090510 at different energy bands detected by the 
Fermi satellite. Panel (b) is the relation of events and time in the 8 - 260 keV energy band, and 
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Figure 1: The light curves of GRB 090510 in different energy bands (figure reproduced from reference Q). 
Panel(a): LAT photons passing the off-line (red) and onboard (blue) event selections. The lines represent 
the relation between photon energy and arrival time for linear (solid lines) and quadratic (dashed lines) 
LIV. Panel(b-f): GBM and LAT light curves, from lowest to highest energies. Panel(f) shows the isolated 
photons with energy F > 1 GeV. 
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the maximum number of events is up to 17,000 counts per second. In Panel (f) there are only 18 
photons with E > \ GeV detected by Fermi-LAT. 

We analyze GRBs detected by the Burst Alert Telescope (BAT) onboard the Swift satellite, and 
concentrate on the 90% of photons in an energy band instead of single photons. BAT can catch 
photons in the energy band of 15 — 150 keV, and record the durations of GRBs, Tqq. It is certain 
that the LIV-induced time delay (Atuiv) cannot be longer than the duration of a GRB. Otherwise, 
the LIV effect would spread the duration of the GRB until it is longer than Afuiv- Therefore, Tgg 
is the upper limit of the LIV-induced time delay. Thus we always have 


Afniv < Tqo. 


(4) 


Eq. dH) is a much more conservative but safer estimation of Atuiv compared to previous works 
using spectral lag as the upper limit. Substituting Eq. Q into Eq. ([3]), we finally get 

(1 -h n)Fn{z) 


Mqg > AE 


2HnE 


0-L90 


where 


Pn{z) = / dz'^ 


(1 + V) 


f\n 


(5) 


( 6 ) 


(1 + z')^ -|- Ha 

From Eq. ([5]), we can see that a rigorous limit on Mqg requires GRBs of short duration. We choose 
GRBs from Swift data archive, and only short GRBs with duration Ego < 2 s are selected. Our 
sample consists of 20 short GRBs with measured redshift in the range z E [0.093, 2.609]. The GRB 
sample is listed in Table [H We use the duration Tgg as the upper limit of Atuiv between lower 
energy (15 keV) photons and higher energy (150 keV) photons. Namely, we get a conservative 
lower limit of the quantum gravity energy scale Mqq. 

Here we mainly focus on the linear correction case, i.e., the n = 1 case. The results are listed 
in the fourth column of Table HI and we get a conservative limit Mqg > 5.05 x 10^^ GeV from 
GRB 140622A for linear correction of the dispersion relation. Most of the GRBs give an order of 
10^^ GeV, and the most rigorous constraint is given by GRB 140622A because of its high redshift 
and short duration. Although GRB 090426 has the highest redshift, its duration is longer than 
most other GRBs. Due to the low redshift and relatively long duration, GRB 061201 gives a limit 
of Mq g on ly at t. 


results |12|.I18 h20|. 


le order of ~ 10^^ GeV. Our constraint on Mqg is much looser than previous 
23-^, because the energy band we used is narrow. However, our result is much 


more statistically significant. This is because the whole energy band in 15 — 150 keV is used in our 
calculation, while other works only chose a limited number of high-energy photons. 
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Table I: The properties of 20 short Swift GRBs we used in the calculation and the resulting quantum gravity 
energy scales, z is the redshift, Tgo is the duration, and Mqg is the lower limit of quantum gravity energy 
scale in the linearly corrected case. 


GRB 

Z 

Tgo 

(s) 

Mqg 

(lO^^GeV) 

150120A 

0.46 

1.2 

0.26 

150101B 

0.093 

0.018 

3.25 

141212A 

0.596 

0.3 

1.34 

140903A 

0.351 

0.3 

0.77 

140622A 

0.959 

0.13 

5.05 

131004A 

0.717 

1.54 

0.32 

130603B 

0.356 

0.18 

1.30 

101219A 

0.718 

0.6 

0.81 

100724A 

1.288 

1.4 

0.63 

090510 

0.903 

0.3 

2.06 

090426 

2.609 

1.2 

1.42 

071227 

0.383 

1.8 

0.14 

070724A 

0.457 

0.4 

0.76 

070429B 

0.904 

0.47 

1.32 

061217 

0.827 

0.21 

2.69 

061201 

0.111 

0.76 

0.09 

060502B 

0.287 

0.131 

1.43 

051221A 

0.547 

1.4 

0.26 

050813 

1.8 

0.45 

2.70 

050509B 

0.225 

0.073 

1.99 


Finally, we also use GRB 140622A to constrain Mqg in the quadratic and cubic correction cases 
of the dispersion relation, i.e., the n = 2 and n = 3 cases in Eq.(l2|). We get a conservative limit 
Mqg > 3.90 X 10^ GeV for the quadratic correction, and Mqg > 3.48 x 10^ GeV for the cubic 
correction. The constraint on the LIV effect in the quadratic and cubic cases is much looser than 
in the linear case. 


IV. SUMMARY 

In summary, we have used 20 short GRBs with redshift measurements from the Swift data 
archive to constrain the possible LIV effect that was predicted by some quantum gravity theories. 






Based on the fact that the LIV-induced time delay should not be longer than the duration of a GRB, 
we derived the lower limit of quantum gravity energy scale at the order of Mqg GeV. The 

strictest limit was given by GRB 140622A, i.e., Mqg > 5.05 x 10^^ GeV. Our constraint on Mqg, 
although not as tight as previous results, is safer and more reliable than previous studies. GRBs 
with higher redshift and shorter duration give a tighter constraint on Mqg. Future observations 
at higher and broader energy bands will also tighten the constraint. 
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